function [f] = fobj(X, Z, V, alphav, sigmax, sigmaz, sigmav)
% % -log(prior)
% f = (0.5/sigmav^2)*norm(V, 'fro')^2 + numel(V)*log(sigmav) + ...
%     (0.5/sigmaz^2)*norm(Z, 'fro')^2 + nnz(Z)*log(sigmaz) - ...
%     logpibp(Z, alphav) + 0.5*(numel(V)+nnz(Z))*log(2*pi);
% 
% % -log(likelihood)
% f = f + 0.5/sigmax^2*norm(X-Z*V', 'fro')^2;

f = -prop(X, Z, V, alphav, sigmax, sigmaz, sigmav);

end